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Abstract 

It is explicitly shown that either the approximate solution of the integral equation 
for the inverse of the pion form facto, r or the result of the Pade approximant method 
of resumming the one loop Chiral Perturbation Theory (CPTH) are equivalent to the 
standard vector meson dominance (VMD) models, using the vector meson coupling to 
two pseudoscalars given by the KSRF relation. Inconsistencies between the one loop 
CPTH and its unitarised version (or the VMD model) are pointed out. The situation 
is better for the CPTH calculation of the scalar form factor and the related S-wave tttt 
scattering. The branching ratios of r — > 7i + tt°u, t — ► Kirv, r —>■ K + 7]v and r — > K + K°u 
using only two inputs as the p and K* masses, or the two corresponding rms radii, agree 
with the experimental data. Using the same number of parameters, the corresponding 
one loop CPTH calculation cannot explain the r data. 
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1 Introduction 



The hypothesis of Vector Meson Dominance (VMD) [|TJ has proved to be an useful and 
convenient concept in low energy hadronic physics. It enables us to describe many low 
energy phenomena, below 1 GeV scale, in a compact and convenient language, although 
not always correct. An apparently different method was based on the dispersion relation 
approach which was the main activities of the soft hadronic physics in the fifties. It was 
soon realised to some authors that the VMD model can conveniently be used to describe 
the more complicated dispersive approach. 

A more recent approach to these problems with a different goal was based on the Chiral 
Perturbation Theory (CPTH) pioneered by Li and Pagels 0, Lehmann ||, Weinberg 
and Gasser and Leutwyler |5], [| [7]. In the Lehmann's approach, the unitarity relation 
of the S-matrix was taken into account but not the chiral symmetry breaking effect. In 
the Weinberg approach, the question of unitarity was ignored giving place to the more 
systematic approach of the perturbation theory in which the unitarity relation was only 
satisfied perturbatively, i.e. order by order, but the chiral symmetry breaking effect 
was taken into account. The latter approach enables us to derive systematic low energy 
theorems, provided that strong interaction can be treated perturbatively, which is an 
assumption and has to be demonstrated. 

In the late fifties, the problem of the validity of the perturbation theory for strong 
processes was questioned, but this issue has been ignored in recent numerous studies of 
Chiral Symmetry using CPTH, either by assuming that the strong interaction involved 
was not sufficiently strong to invalidate the perturbative approach, or that the CPTH 
was effectively a low energy power series expansion in momentum which could evade the 
unitarity constraint. 

There are few publications on Chiral Symmetry where chiral symmetry breaking ef- 
fect and unitarity are simultaneously taken into account. This last approach, combining 
the CPTH with Unitarity, will be called as the Unitarised Chiral Perturbation Theory 
(UCPTH). It consists in deriving a similar expression as the one loop CPTH result but is 
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supplemented by using the inverse amplitude, the N/D or the Pade approximant methods 
in order to satisfy the elastic unitarity relation [|3], |8|, |9], [H], II]. This procedure enables 
us to extend the perturbation theory to incorporate low energy resonance or bound state 
phenomena which the standard CPTH cannot handle. We can, in principle, include the 
inelastic effects but calculation becomes less simple. Except for a most simple calculation 
of the inelastic effect, we deal in this article the low energy phenomena where the elastic 
unitarity relation dominates the calculation. 

The elastic unitarity relation for the form factor, which is incorporated in the UCPTH 
approach, leads to the following two consequences. First, the expression for the form 



factor is an integral equation of the Muskhelishvilli-Omnes type Jl2| |T3[ which should 
not be solved perturbatively. Second, the phase of the form factor in the low energy 
region (where the inelastic effect can be neglected), must be the same as that given by 
the corresponding strong interaction phase shift J2TJ. The solution of the Muskhelishvilli- 
Omnes integral equation |12| has a polynomial ambiguity, to be discussed below, but the 
phase theorem is independent of such an ambiguity. Comparision between theory and 
experiment should therefore be done with the magnitude and phase of the form factor. 

The one loop CPTH calculation for the vector and scalar pion form factor was done 
a long time ago ||. A 2 loop CPTH calculation was also recently done |14|, These 
papers treated the form factors perturbatively and therefore cannot be used to calculate 
the r decay to 2 pseudoscalars which is dominated by the p and K* resonances. 

In a recent paper on the study of the isovector form factors ||, it was shown that the 
question of the unitarity must be respected in order to describe in a simple manner the 
low energy pion physics with or without resonances. It was explicitly pointed out that 
two approaches could be used: a) One could use the standard CPTH but after having 
done the one loop perturbative calculation, one must resum the series by the Pade method 
in order to satisfy the elastic unitarity relation []3], |], |9|, |H], |Tl| (this approach could be 
regarded as the large number of flavor Nf expansion). Other applications of the Pade 
method in physics have been shown to be successful [|T(J, but we give here the reason 
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why it should be. b) One could write the dispersion relation of the inverse of the form 
factors or the scattering amplitudes which takes automatically into account of the elastic 
unitarity relation, and then solve approximately the resulting integral equation || [T7| . 



The well-known N/D method could be regarded as an extension of this method but is 
more flexible. 

Both methods yield the same result and give rise to the well-known KSRF relation for 



the width of the vector meson |n| . It was implicitly shown that these two methods are 
equivalent to the VMD model || |T|, f|J using the pirn coupling as given by the KSRF 
relation ||18|| . The present method contains less parameter than the VMD model because 
the KSRF relation is a direct consequence of our approach and not an added assumption. 

In this paper we want to point out the equivalence between the VMD model on one 
hand, and the Pade, the inverse amplitude and the N/D methods on the other hand, 
for the calculation of tctt, KK, ttK and rjK vector form factors. They have the same 
number of parameters as those used in the one loop CPTH method but yield a much more 
significantly improved prediction due to the respect of the constraint of the unitarity (by 



a factor of more than 50 in the nK form factors squared) |8|, [L9|, |20 



We want to improve the inverse amplitude method, by expressing the form factor in 
terms of the N-function of the N/D method for the corresponding strong amplitude. We 
briefly give the reason why the N function used in the vector pion form factor calculation 
can be approximated by the Weinberg low energy expression. 

As mentioned above, the unitarity relation for the form factor would require that its 
phase must be that of the corresponding strong interaction, if the inelastic effect was 
unimportant in the physical region of interest |21| . This is the only model independent 
relation which we can make, apart from the (more fundamental and rigourous) Ward 
identity on the the non renormalisation of the charge of the pion at the zero momentum 
transfer. For the vector and scalar pion form factors, because there are no significant 
inelastic effect below 0.9 GeV 2 , the phase theorem should be valid in this region. Hence 
the phase of the vector pion form factor should have a phase of 90° at 770 MeV, the p 
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vector meson mass, and similarly, the 7rK vector form factor has a phase of 90° at the 
K* mass. Throughout this article and in the UCPTH approach, we have to rely on this 
phase theorem to check the validity of our approximation scheme. This point is ignored 
in the standard CPTH calculations. 

As a direct application of the UCPTH method, we show that the branching ratios of 
r — > Vector Mesons + v agree reasonably well with the experimental data using only 
inputs as the p and K* masses, or alternatively, the corresponding rms radii ( which the 
usual CPTH cannot be used or is wrong, being too low, by a factor of 10 to 50). The 
branching ratios of r — > K + i]u and r — > K + K°u are however smaller than the observed 
branching ratio by a factor of 2. A much better agreement with the data is obtained 
by introducing the polynomial ambiguity as a simulation of the inelastic effect, e.g. the 
contribution of the urn intermediate state in the unitarity relation. By fitting the height 
of the pion form factor data at the p resonance with the first order polynomial, the r.m.s. 
radius is now in good agreement with the data. Alternatively, using the experimental 
r.m.s. radius as an input, the magnitude of the pion form factor at the p resonance is in 
a very good agreement with the experimental data. 

We also want to show that the inelastic effect of the KK states on the pion form factor, 
the inelastic effect of the r/K intermediate state on the nK form factors are however 
negligible. A more correct approach, using the coupled two-channel unitarity calculation, 
will be done in the future. 

From a more fundamental viewpoint, we question the assumption made in the CPTH 
approach that perturbation theory can be used to study the vector pion form factor which 
could involve important non-perturbative effects. For this purpose, we want to point out 
some novel theoretical inconsistencies between the one loop CPTH and its unitarised 
version UCPTH calculation of the vector pion form factor (or equivalently the VMD 
model); the same observation is also valid for the P-wave im elastic scattering. Our 
observation on the inconsistency can be extended to the two-loop CPTH calculation, and 
is conjectured to higher orders. We point out, however, that the problem of unitarity may 
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not be serious in the CPTH calculation of the low energy 1=0 scalar form factor or the 
related elastic S-wave tttc scattering. 

It is our feeling that, not only because the Breit Wigner form cannot be expanded in a 
convergent power series of momenta near the p resonance, because of this inconsistency, 
even with higher order loop calculation, CPTH will remain as a incomplete low energy 
perturbative theory and cannot handle resonance or bound state problems which are a 
manifestation of the non perturbative effects. 

The plan of this paper is organised as follows: In section 2, which is the main part of 
our paper, we give a detailed calculation of the vector pion form factor by perturbation 
method and also by the non perturbative inverse amplitude and the N/D methods where 
the elastic unitarity is satisfied. In section 3 the scalar form factor is calculated also by 
two methods. We then point out that the use of the perturbation method can be justified 
here. In section 4, the equivalence between the UCPTH and the VMD is shown. Section 
5 deals with the SU(3) generalisation of the form factors. Section 6 is devoted to the 
applications of the above calculations to the r decays. 

2 Calculation of the Vector Pion Form Factor 

a) Theoretical Consideration: One Loop CPTH Result 

We begin first by recalling the main features of the vector pion form factor calculation 
using CPTH || ||. Let us define the vector pion form factor as: 

< n + ( Pl )n-(p 2 ) | V*(0) | >= i{p2-pi)»V{s) (1) 

where s is the momentum transfer squared s = (pi + p<i) 2 and ^(0) = 1. Using the 
dimensional regularisation scheme, the one loop CPTH result is [Q: 

V^{s) = 1 + 2 -^s + ^L_((s - 4m>)H m (s) - sLog^ - |) (2) 
where s is the the momentum transfer, ag(/i) is one of the renormalized constants 



defined by Gasser and Leutwyler 0, ^=93 MeV, the experimental value of the pion 
decay constant, and 



Is - Ami y/s + Js- 4m 2 n U _ 4 m 2 

tf mW = (2-2^_ -^Log £- (3) 

for s > Am 2 ,; for other values of s, H^^s) can be obtained by analytic continuation. 
We can explicitly introduce the expression for r.m.s radius of the vector form factor in 
Eq. (H). Using the definition ^'(0) = \r\r = l/s Rl we have: 

1 2«g(/i) 1 , T m 2 

7ir^r-^ii (Lo9 v +1) (4) 

Using this expression in Eq. (g), we have: 

V pert is) = 1 + — + wAr 2 ((s - ^ml)H^(s) + ^) (5) 

It should be noticed that the real part of the sum of the last two terms on the R.H.S 
of Eq. (H) behaves at low energy as s 2 . In using the r.m.s. radius as an experimental 
input to renormalise the calculation of V pert '(s), we can calculate perturbatively the vector 
form factor to the order of s 2 (modulo of a logarithm). In the language of the dispersion 
relation approach, the r.m.s radius of the pion vector form factor is used as the subtraction 
constant. 

b) Non Perturbative Approach: Elastic Unitarity Constraint 

At issue here is whether the perturbative hypothesis used in the CPTH approach 
justified. We shall discuss this problem later. At the moment, we shall discuss the 
non perturbative problem by using the integral equation approach. Let us review some 
fundamental properties of the form factor V(s). Because of the Ward identity and the 
analytic property of the form factor, we can write the following disperssion relation for 
V(s): 

V(s) = 1 + — + - r a{z) . ^ dz (6) 

S Rl 7T J Ami z \ z ~ s — lt ) 
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where an extra subtraction is made for convenience, and the spectral function o~(z) 
must be taken to be real from time reversal invariance. 

The Muskhelishvilli-Omnes integral equation \Y1\ is obtained by using the elastic uni- 
tarity condition for the spectral function: 



v(s) = i + j- + 1 r y(^>s t ,Mz) dz 

s Rl tc Jiml z I (z — s — ie) 
where 8\ is the P wave tttt phase shift and we have made 2 subtractions in the dispersion 
relation for later purpose. Because of the of the reality condition of the spectral function, 
in the region where the elastic unitarity relation was valid, the form factor would have 
the same phase as that of the strong amplitude. This is the content of the Watson phase 
theorem [^1] . 



The solution for the integral equation (0) is well-known [|I^] and can also be obtained 
perturbatively by an infinite iteration of the integral equation: 

V(s) = P n (s)exp(- r 6l(z) dz) (8) 

7T J Am* Z \ Z — S — 16) 

where P n (s) is the polynomial ambiguity with real coefficients, normalising to unity at 
s=0; they could represent, at low energy, uncalculable higher energy inelastic effects. We 
shall assume at the moment that P n (s) = 1 in the physical region of interest. As can be 
seen, the phase theorem due to the unitarity of the S-matrix is automatically satisfied. 
An extension of this integral equation to include a phenomenological inelastic spectral 
function was previously considered and its solution is also known [|13[ (see below). 

The solution of the integral equation could have been guessed from the reality condition 
of spectral function: As mentioned above, from this condition, the phase of the form factor 
is the phase of the strong amplitude. Because e lS is not an analytic function, the only 
correct solution for V(s) is given by Eq. (^). The one loop CPTH solution, Eq. 
is the once iterated solution of the integral equation Eq. (0): This perturbative result 
can be obtained by setting V(s) in the integral of Eq. (|5|) to be unity and f( ree (s) = 
e tSl ^ sin5i(s) / p(s) — (s — Am^) / (96% f%) which is purely real; this result follows from the 
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evaluation of the one loop Feynman graph using the Cutkosky rule for the absorptive part 
and dispersion relation. The once iterated solution of the Muskhelishvilli-Omnes integral 
equation makes no sense, because the exact solution of this integral equation, as tated 
above, requires an infinite iteration. This problem was discussed in reference |§. 

In order to compute the pion form factor, we can either use the experimental P-wave 
phase shift 5\ or the calculated 8\ phase shifts (obtained from the construction of a 
unitarised P-wave amplitude) in Eq. (H). To see the effect of the exact solution of the 
integral equation, let us observe that the P-wave irn phase shifts passing through 90° at 
the p mass, hence we approximate the P-wave phase shifts 5 by tt6(s — s p ), the zero width 
approximation for the p resonance; the exact solution obtained from Eq.(H) is therefore 
Sp/ (s p — s). It is clear from this example that we should develop a power series expansion 
for the inverse of the form factor, if we wanted to incorporate the low energy property of 
CPTH and at the same time a non perturbative result to take into account of a possible 
resonant behavior. This was exactly done in reference || because the inverse amplitude 
has a nice property that the elastic unitarity relation for the tttt — > tttt partial wave is 
automatically satisfied. 

The non-relativistic limit of this expansion was done a long time ago and is known 



as the "effective range theory" p3| ; it is a momentum power series expansion of kcot5 
which preserves automatically the unitarity of the S-matrix. This is so because the elastic 
unitarity of the partial wave enables us to write f(k) = e lS sinS/k = l/k(l — icotS) , 
where k is the c.o.m. momentum, hence the power series expansion in kcotb which is 
the effective range expansion, is an expansion of the inverse of the amplitude. This type 
of expansion enables us to handle very well the low energy bound state (triplet S-wave) 
and the resonance (singlet S-wave) of the low energy nucleon-nucleon scattering. The 
missing part of the Hilbert transform or analyticity is presumably unimportant for a 
non-relativistic theory. ( Unfortunately, standard explanations of the low energy effective 
range theory do not emphasize the question of unitarity of the S-matrix which is much 
more important than the potential shape dependence of the expansion). 
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As applied to the form factor problem, to the extent that the discontinuity of the left 
hand cut for the partial wave amplitude can be neglected, or can be treated perturbatively 
(see below), the pion form factor can straightforwardly be obtained by considering the 
integral equation for the inverse form factor amplitude ||. This result is equivalent to 
resum the perturbation series using the infinite bubble summation of the pion loops. 

Let us now carry out our analysis, without making any assumption on the left hand 
cut structure. Because the form factor has a cut from 4m^ to infinity, so does its inverse, 
apart from the contribution coming from the zeros of the form factor appearing as poles 
which we assume to be absent here. We shall later give a phenomenological description 
of these zeros as uncalculable inelastic effect which has not been taken into account here. 

The inverse of the form factor V _1 (s) satisfies also the same analytic property as the 
form factor V(s), hence we can write a dispersion relation for it. Assuming that V(s) 
does not have a zero or its position is far from the physical region of interest, using the 
elastic unitarity condition, we have: 



V-\s) = l / y —^dz (9) 

Sj?! VT J Ami Z A \Z — S — It) 

From the general property of the analytic function, the partial wave amplitude f\{s) = 
e l5l sinSi/ p(s) where p(s) = yl — 4m^/s, fi(s) can always be written as the product 



00 V*- 1 {z)e~ iSl ^sinS 1 (z) 



of the 2 cuts, the right hand cut or the unitarity cut and the left hand cut due to the 
exchanged contributions: fi(s) = N(s)/D(s). Because we use the hypothesis of the elastic 
unitarity, the right hand cut function D~ l (s) has the same phase representation as given 
by Eq.(||) with the normalisation D _1 (0) = 1 and where we set P n (s) = 1. Hence we can 
set V*~ 1 (s)e~' lS ^ s 'sin5(s) = p(s)N(s). Eq. @ can now be rewritten as: 

v-\s) = 1 - j- - s - r do) 

S Rl 7T Hml z z {z — s — ie) 

This equation expresses the pion vector form factor in terms of the N function of the 
elastic P-wave irn scattering amplitude instead of the P-wave phase shift 5i as given by Eq. 
(^). For the following purpose of calculation, Eq. ( |TD| ) is more convenient. As mentioned 
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above the function N(s) involves a dispersion integral of the discontinuity of the partial 
wave amplitude across the left hand cut involving the effects of the / = 0, 1 and 2 pair 
of pions exchanged in the crossed channels. Because of the normalisation of the function 
D(0) = 1 and that the Weinberg low energy theorem is operative at low energy for the 
elastic tttt scattering amplitude, we can write the N function for the P-wave as: 

*M = i£ H^#{ i + sftrf^^ r , ': N 2 \;;\ (id 

96nf^ 7r Jo [z + Aml) z {z + s) 

Eq. ( |10D and this equation are still exact, apart from the assumption on the absence of 

the zeros. It provides an alternative solution to that given by Eq. (H) on how to construct 
the solution of the form factor problem when the dynamics of the strong interaction is 
known. We now discuss an approximate scheme for Eq. (0). To calculate the form factor 
V(s) in Eq. ([[(]) we must know the function N(s) for s > 4m^. In this energy range, the 
contribution from the ImN(s) is usually small because the denominator in Eq. ( |TT| ) never 
vanishes. Hence we approximate N(s) by the first term on the R.H.S. of Eq. (JTTj) , i. e. 
we neglect the contribution from ImN(s) or equivalently to use the Weinberg low energy 
theorem for tttt scattering [^]. The approximation scheme is therefore to start initially 



with the Weinberg low energy theorem for the partial wave amplitude, then this result 
can be improved by correcting for the deviation of the Weinberg theorem at higher energy. 
This procedure is reasonable because we use a twice subtracted dispersion relation in s 
which emphasise the low energy contribution at small value of s and hence it is valid to 
approximate the nir amplitude by the Weinberg expansion as a first approximation. 
The form factor V(s) can now be written as: 

V{S) = 1 - s/s Rl - ^{{a - AmDH^is) + 2s/3} (12) 
Eq. (|12[) can also be derived by using the Pade approximant method: 

T/free 

= y iooP (13) 

ytree 



11 



where V tree refers to the tree amplitude which is equal to unity and y x - lo °v refers to 
the one loop amplitude, i.e. the last two terms on the R.H.S. of Eq. (^). The Pade 
approximant method yields the same result as that given by Eq.fll^D which satisfies the 
elastic unitarity relation. 

We can justify the use of the Pade method to resum the one loop perturbation series 
by looking at the imaginary part of the inverse of Eq. fll3l): it is equal to —p(z)N(z) 
with N(z) given by the tree amplitude and hence containing no left hand cut, or N(z) = 
(s-4ml)/(96nm. 

The problem of tttt scattering, in the approximation where the contribution of the 



left hand cut could be neglected, was done a long time ago by Brown and Goble [25 



The more exact calculation with the contribution of the left hand cut and in the chiral 
limit was done by Lehmann [|3| where the real part of the logarithm terms due to the 
right and left hand cuts cancel each other out in the chiral limit. A more exact one loop 
perturbation calculation was done by and its unitarised form was given by |]] . In this 
last work, the calculated P-wave phase shifts agree very well with the experimental data 
from the inr threshold to lGeV. and differs little from those calculated by neglecting the 
left hand cut contribution. 

Similarly, our calculation for the form factor using Eq. ([12]) differs by at most 1-2% 
from the direct calculation of the form factor using Eq. (|||), with the P-wave phase shifts 
given by the unitarised version of the P-wave irn scattering || . Explicit results are given 



in reference ||22|| . Hence to the extent that the pion loops are taken into account, the left 
hand cut contribution to the N-function of the P-wave amplitude can be neglected. 

From another point of view, the problem of treating the exchange of p and a as particles 
with their full propagators, in the framework of the generalised linear sigma model, was 
previously examined in the reference p6 |. It was found that, in order to have the validity 
of the KSRF relation, and treating the problem at the tree level approximation, the effect 
of p and a exchanged in the t, u channels should cancel out approximately, when their 
masses satisfy the relation m a = \/2m p which is not in contradiction with the experimental 
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fact. This problem should be further investigated. 

The above discussions give a justification for the approximation done by Brown and 
Goble [25|. In fact Eq. flT2p is exactly the inverse of the D-function which is the same as 



the Omnes function given by Eq. @. 

In this section, we limit ourself to the approximation where the discontinuity of the left 
hand cut can be neglected in order to make a connection of our calculation scheme with 
the perturbation series and also with the large number of flavor Nf expansion scheme. 
Objections can be raised against this approximation because of the large violation of the 
crossing symmetry which we were taught back 30 years ago to be an important property 
of the field theory. This opinion is certainly correct, but calculations done at that period 
were without the constraints of chiral symmetry, i.e. without having the low energy 
theorems which can be used to write subtracted dispersion relations. Subtracted dispersion 
relations are important because we want to suppress the high energy contributions which 
are difficult to calculate. With the use of the low energy chiral theorems in dispersion 
relations, the crossing symmetry difficulties are minimized as the scale of the physical 
problem at low energy is fixed by the low energy theorems. 

c) Deviation from the Elastic Unitarity Constraint 

As we shall show below, the inelastic contribution in the unitarity relation for the 
vector pion form factor due to the Kaon loops amounts to a few percent. The situation is 
quite different when we consider the contribution of the ujtt and possibly the pirn channels 
which are experimentally important. The calculation of this contribution to the pion form 
factor is quite complicated. We can derive a similar equation to the Muskelishvilli-Omnes 



equation |13[ but the phenomenological application has considerable uncertainties |27 . 
Instead of the integral equation (|?p, we have now: 

v W = 1 + -L + ^r VMW + 'M dz (14) 

S Rl TT J 4ml Z z {Z — S — lt) 

where fi(s) = — tj being the inelastic factor, and crj(s) being the inelastic contri- 

bution to the unitarity relation and differing from zero for s above the inelastic threshold 
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Si. The solution for this equation is 

We also have to require that the first derivative of V(s) at s=0 to be sr-i. 

Below and sufficiently far from Sj, e.g. the ujtt threshold, we can roughly parametrise 
the contribution of the second term on the R.H.S. of Eq. ([y|) as a polynomial P n (s) which 
is real for s < Si and is the polynomial ambiguity in the solution of the Muskhelishvilli- 
Omnes equation. Because of the normalisation P n (0) = 1, the introduction of this factor 
will not influence the Ward identity, but does influence the value of the r.m.s. value of 
the pion radius. This type of approximation is reliable for energy below the inelastic 
threshold, but is erronous at the inelastic threshold and also at higher energy. 

We shall fit the experimental data below with the simple expression (l + as/s p ), where 
a is a constant in order to simulate phenomenologically the inelastic effect. The pion 
form factor can now be written as: 

V(s) = 1 + as / s P (16) 

^ ' 1 - s/s Rl - ^{(s - 4m*)H m (s) + 2s/3} 1 ; 

d) Phenomenological Applications 

1) Let us first calculate the form factor without using the phenomenological introduc- 
tion of the inelastic effect in the unitarity relation. Eq. fll^D gives then the expression for 
the vector pion form factor. Because the vector r.m.s. radius is positive, the vector form 



factor has a resonant character. Its width satisfies the KSRF relation ||18|| . The relation 
between the p resonant mass squared s p and is: 

SRl = 1 - 96^{(sp - 4ml)ReH w „(s p ) + 2s p /3} 
In Eq . (|i~7l) , we can either use s Rl or s p as an input parameter. If the r.m.s. radius 
was used then one would predict the p mass and width which can be seen to satisfy the 
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KSRF relation [18|: T p = (s p — 4m 2 r ) 3 / 2 /(967r/' 2 ) where we have neglected the finite width 
correction [p9 |. In terms of the g p7T7T coupling constant, we have: 
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2 S P 



where g pn7T is defined at s = s p . Using the experimental value m p = 0.77GeV, the 
uncorrected width of the vector meson p is 0.141 GeV. Including the finite width correction 
|2T| , the p-width, defined as {^fSpT p )~ l = cot5i(s)' \ s=Sp , where the sign ' denotes the first 
derivative of cotSi with respect to s. The KSRF relation, with the finite width correction, 
yields T p = 155MeV which is very near to the experimental value of 151.5 ± 1.5MeV. 

If one used s p as an input parameter, then one would predict the p width by the 



KSRF relation and the r.m.s. radius of the pion. For a number of reasons which 
will become clear later, we use the p mass as an input. The calculated r.m.s. radius is < 
r v >= 0.40 ± 0.01/m 2 compared with the experimental value < r v >= 0.439 ± 0.008/m 2 



| 30fl . The calculated value is therefore too low. This difficulty is related to the fact that 



the calculated pion form factor at the p peak is also too low, as can be seen in Fig. (1) 
(see below). This result is expected, because we neglected the inelastic contributions in 
the unitarity relation which does not change the imaginary part of the pion form factor 
but does contribute to its real part as explained above. A much better agreement with 
the experimental data, both at the p peak and also for the pion r.m.s. radius can be 
obtained by phenomenologically introducing the inelastic effect. 

The phase and modulus squared of the pion form factor calculated by the CPTH, Eq. 
(|), and by UCPTH, Eqs.([]]|[l3|), using < r\ >= 0.40 ± 0.01/m 2 are shown on Figs. 
(1) and (2). Although they both agree with each other at low very momentum transfer, 
due to the dominance of the r.m.s radius term, the higher energy experiment data are 
much better represented by the UCPTH calculation than those of the CPTH calculation. 
Experiences with analytic functions show that a small difference between two analytic 
functions in one region can be greatly amplified in another region. This explains why the 
form factors calculated in both CPTH and UCPTH have the same value and radius at 
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s=0, but have completely different behavior around the p resonant region. The UCPTH 
calculation has the p resonance character while that from the one loop CPTH does not. 
The phase theorem of the pion form factor PTJ| can be used to test whether the unitarity 



relation is satisfied or not. From Fig.( 2), for s > 0.2GeV 2 the CPTH phase is clearly 
much less than those obtained experimental data reflecting that the unitarity relation is 
badly violated in this approach. The UCPTH pion form factor phase is, on the other hand, 
in very good agreement with the data from the 27r threshold to lGeV 2 . This reflects that 
the UCPTH approach satisfies the unitarity relation in this region. We conclude that the 
unitarity relation plays a crucial role in the low energy calculation of the pion form factor 
and also of the tttt scattering. We shall show elsewhere some difficulties of calculating the 
process ir — > 77* encountered in the CPTH approach can also be removed by using the 



UCPTH method 123 



e) Inadequacy of CPTH calculation for Vector Form Factor. 

Let us now examine Eq. ([5]). As noted above, the last term in Eq. (|5]) has a different 
low s behavior for the real and imaginary parts: at small s, its real part behaves as s 2 
while its imaginary parts, in the chiral limit, behaves as s. This mismatch between the 
small s behavior of the real and imaginary parts is due to the perturbative approach which 
results in a violation of the phase theorem. In order to restore this theorem, even at low 
energy, the next order term in the imaginary part has to be included in the perturbation 
calculation as pointed out in ref. || and will be discussed later in the scalar form factor 
calculation. This problem does not exist in UCPTH approach as can be seen from Eq. 



(0) 



From the discussion of the previous paragraph, we are led to the study of the 0(s 2 ) 
term in Eq. ([5]) and that in Eq. ([121) . Although they are numerically very small at 
low values of s, because we are interested in extending the region of the validity of the 
perturbation theory to near the p resonant region, they must be taken into account in 
this energy range. Hence we are led to examine the consistency of the CPTH approach 
compared with the UCPTH calculation or the VMD model which we shall show later to 
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be equivalent. For this purpose, not only we want to compare the experimental data on 
the modulus of the pion form factor but we also want to use the phase theorem which 
is model independent, to compare the calculated phases with those obtained from the 
experimental data. 

From Eq. fll2|) , it is seen that the perturbative approach of Eq. (||) leaves out, at low 
energy, a term {s/s^) 2 which should be compared with the real part of the last term in 
Eq. (||) because they are both of the same order in s. The perturbative approach is only 
justified if the last term is much larger than the former. Let us now define the parameter 
ry, defined as the ratio of the real part of the last term on the R.H.S. of Eq. @ over 

(*/^) 2 : 



((s - AmDReH^s) + f) 



w'Lr " (19) 



In Fig. (3), the ratio ry is plotted as a function of s. It is seen that, over a wide range 
of values of s, including small values of s , ry is much less than unity instead of being 
much larger than unity in order to guarantee the validity of the perturbation theory. As 
we shall show below the situation is very different in the calculation of the scalar form 
factor. 

We have discussed so far, the CPTH one-loop calculation of the pion form factor. It is 
not difficult to calculate the two-loop contribution to the vector pion form factor as was 
previously done numerically by Gaisser and Mei/3ner |L4] and more recently analytically 



by Colangelo et al. |jTq| . The main point is that the two loop calculation requires an extra 
subtraction in the dispersion relation and hence it is necessary to introduce an extra 
parameter, the second derivative of the pion form factor at the origine. These authors fix 
this parameter by taking the VMD model which is (s/s p ) 2 . What one calculates in the 
CPTH approach is the contribution of the s 3 terms, modulo of the logarithm. One can 
then show the result of the calculated 0(s 3 ) term is still smaller than the {s/s p ) 3 terms 
coming from the VMD model, although the corresponding ry term is, on the average, a 
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factor of 2-3 larger than the value of ry of the one loop-calculation. 

Although we have not carried out the calculation of the 3-loop and higher loop am- 
plitudes, we suspect the same difficulty also occurs. This unsatisfactory situation of the 
1-loop, 2-loop and possibly higher loop calculations, is due to the use of the perturbation 
theory to study nonperturbative phenomena. One would have to increase the p mass to 
be larger than 1.5 GeV in order to make accuracy of the vector form factor calculation 
comparable to the scalar case(see below). 

It is so much simpler and more transparent to use Eqs. (0,0) which is even better 
than the VMD model as a starting point. One then improves it by taking into account 
of the 2-loop graphs which are neglected in the calculation. A better expression, which 
takes into account of the inelastic effect yielding a correct r.m.s. radius and the height of 
the p peak Eq. is even better than Eqs. ( |12] , |T3j ) for this purpose. If the major goal 

of physics was to describe experimental data by a few parameters, because CPTH had 
to rely on the VMD model to get the necessary counter terms, it would become a rather 
special way of doing physics. 

We would like to ask whether by including higher order terms in the CPTH approach 
that one could generate the Breit-Wigner form. We have a strong doubt about this 
possibility. Our feeling is based on the experimental fact that there is very little inelastic 
effect below 1 GeV in the form factor so that higher loop effect cannot be important in 
the p resonance region. 

In Fig. (4) we plot the modulus of the pion form factor calculated by the UCPTH 
compared with the sum of the first four terms of the Taylor series expansion of the UCPTH 
calculation. It is seen that, at low energy, the first four terms are a good approximation 
for the experimental data but we cannot make the the Breit Wigner form i.e. to make the 
high energy curve to turn down. This illustrates the problem of the non convergence of 
the perturbation series for the amplitude; there is, however, no difficulty in generating the 
Breit Wigner form by expanding the inverse amplitude as a Taylor's series of momenta. 
Our remark here is a relativistic generalisation of the Bethe's "effective range theory" to 
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handle the resonance and bound state in the two nucleon system |23| . 



2) We want now to improve our calculation by taking into account of the inelastic effect. 
As explained above, the simplest parametrisation can be done by setting P n (s) = 1+as/ s p 
and is given by Eq. fll6|). 

The good fit to the data can be obtained with a = 0.15 with ^/s^ = 0.77GeV. The 
vector pion r.m.s. radius squared is now < ry >= 0.46 ± 0.01/m 2 compared with the 
experimental value < r v >= 0.439 ± .008/m 2 [3C]. A slightly better fit to the data is 
obtained with a = 0.13 but the p mass is now changed to 0.773 GeV; the vector pion 
r.m.s. radius squared is now < r\ >= 0.45 ± 0.01/m 2 which is quite good. 

Before leaving this section let us write down the dispersion relation for the form factor 
when the r.m.s. radius is not used as a subtraction constant, but a low energy measured 
form factor at s = — Sq is used as a subtraction constant. Let us consider a measured 
space-like form factor at s = —sq (the time-like form factor can straightforwardly be 
written down). Instead of Eq. (|7j), we now have: 

V {S ) = V(0) + (V(0) - V(-s ))± + *±±!> f a (20) 

S 7T J4rnl Z{Z + S ){Z — S — ie) 

The one-loop CPTH result can be obtained perturbatively from this equation by setting 
in the integrand V(z) = 1 and fi ee (s) = e tSl ^ sinSi(s) / p(s) = s/(967r/ 2 ). Because we 
use the input as the measured form factor at s = —sq, we cannot discuss the obtained 
results as an expansion in a power series in s. 

Let us finally discuss the ghost problem in the Pade, inverse or the N/D methods. 
Eqs.([T2|,[!3|) develop a pole at s g = -3.8 • lO^eV" 2 which is a ghost and should not be 
there. It has to be eliminated from these equations by multiplying them with the factor 
(1 — s/ s g ) which does not effect the normalisation of the form factor at s — 0, but change 
the value of the pion form factor by an amount s/s g which is completely negligible at 
small s. The phase theorem is also unaffected by this factor. We show here that the 
presence of a ghost due to the unitarisation procedure can be tolerated, as long as it is 
sufficiently far from the physical region of interest. 
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3 Calculation of the Scalar Form Factor 



Similarly to the calculation of the vector pion form factor, we can now calculate the 
isoscalar scalar form factor S(s) defined as: 

< 7i a {p 1 )7T b {p 2 ) | m(uu + dd) I >= ml5 ab S(s) (21) 

where s = (pi + p<i) 2 . After introducing the scalar r.m.s. radius, using the definition: 
S (0) = I < r i >= V s /?2 to eliminate the scale dependent fi, the one-loop perturbative 
result for the scalar form factor is: 

S pert Xs) = 1 + — + jAr 2 ((s - ^/^(s) + ±) (22) 
After unitarisation, the UCPTH result for the scalar form factor is: 

S(S) = -^((s-ml/2)H mr (s) + ^) (23) 

Similarly to the definition of ry, we can define r$, the ratio of the 0(s 2 ) of the CPTH 
calculation and (s/sr 2 ) 2 : 

rs = k n V2 24 

In Fig. (3), the ratio r s is plotted against the energy squared. It is seen that r s 
is larger than unity, although not much larger on the average; hence, we are assured 
about the approximate validity of the perturbation theory. The reason for the difference 
with the vector form factor calculation is due to the coefficient of the chiral logarithm 
term being larger by a factor 6 in the scalar form factor. This fact explains the partial 
success of the one loop CPTH calculation in the S-wave nir scattering ( this was done with 
the prescription that the phase shift is proportional to the real part of the partial wave 
amplitude in tttt elastic scattering ||. Had the calculation been done with the definition 
tancj) = Imf/Ref, the result would be very different). 

It should be remarked that in the CPTH approach, the calculated phase of the form 
factor is not the same as the phase shift S because the unitarity relation is not satisfied 
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in the perturbative approach. In fact, if we calculated the phase of the form factor in the 
CPTH approach, Eq. (|22|), we would find out that: 

p(s)(s — m 2 12) 

tan{<p CPTH ) = 167r/ 2 {1 + S/SR2 + _j_ [{s _ m 2 /2)ReHinr{s) + ( 25 ) 
to be compared with the phase calculated by the UCPTH method, Eq. (|23[) : 

, / , x _ P(s)(s-m 2 j2) 

tan(cf> UCPTH ) - i6nmi s/s ^ __j_ [{s _ ml/2)ReHUs) + X]] ( 26 ) 

In Fig. (5) the dashed curve is the result of the standard one loop CPTH using the same 
parameter as that of UCPTH Eq. (^|). The phases of the scalar form factor calculated 
by CPTH and UCPTH, 4>cpth and 4>cpth are plotted against the energy squared. It is 
seen the phase of the form factor calculated by CPTH is in a much better agreement with 
the data than the vector case; the agreement between the UCPTH calculation and the 
experimental data is, however, better than the phase of the one-loop CPTH calculation. 
It should be recalled that the phase shift calculated by the CPTH method was done by 
making the prescription that it is proportional to the real part of the strong S-wave irn 
partial wave and agrees well with the experimental data. The difference between the 
CPTH phase of the form factor and the phase shift is therefore a measure of the violation 
of the unitarity relation. 

It is seen that these two phases are different even at the threshold. Let us compute 
the limit s — > 4m 2 for tanip/ p(s). From these two equations and use the definition for the 



scattering length a, the limit s — > 4m 2 of tancp = (l/2)ys — 4m 2 a, we get: 

_ 7m„ Ami , x-i (97] 

acPTH - — „ i H h 2 J [2( 

32nf 2 s R2 247r 2 f 2 

and 

Im-x Ami 1 

aucPTH = „ f2 (1 - . 2 n ) (28) 

For sr 2 = 0.56GeV 2 which agrees with the experimental scalar radius < r 2 >= 
0.41/m 2 , we have: a C pTH = 0.123m" 1 and Qucpth = 0.23m" 1 . These values are to 
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be compared with the direct CPTH calculation of tttt scattering a = 0.20m n l H and the 
same process using UCPTH, a = (0.26 ± 0.02)m- 1 §. 

It is clear that the phase theorem for the form factor is violated in the CPTH method. 
How can we rectify this situation for CPTH method? As we mentioned above, there is 
a mismatch between the real and imaginary part of our calculation. In the perturbative 
calculation, the real part of the form factor is calculated to the second order in s while its 
imaginary part is only calculated to the first order in s. In order to compute correctly, in 
the CPTH approach, the phase of the form factor given by tancj) = ImS/ ReS, we must 
also calculate the imaginary part of the form factor to the second order in s, i.e. the two 
loop order for its imaginary part. More precisely, the perturbative unitarity relation for 
ImS including the two loop contribution is: 

ImS(s) = p(s){S 1 (s)f 1 (s)* + S 2 (s)f*(s) + S 1 (s)f;(s)} (29) 

where Si and £2 are, respectively, the first and the remaining terms on the right-hand 
side of Eq. (^); they represent the tree and one loop amplitudes. Similarly fi and / 2 
are the correponding elastic tttt scattering amplitudes. The calculated scattering length 
is now : 

7m w , Refn(Amt) , 4ml llml . . . 

Using the value -Re/2(4m^)//2e/i(4m 2 ) ~ 0.25 from the one loop calculation for the 
1=0 7T7T scattering amplitude || |9j, this Eq. yields a = O.lOSm^ 1 . This value is very 
close to the scattering length a = O^Om^ 1 obtained by CPTH method || for the elastic 
tttt scattering using the prescription that the phase shift 5 is equal to real part of the 
calculated partial wave. 

We see that by correcting for the low energy mismatch of the real and the imaginary 
parts we can get a reasonable agreement with the phase theorem. A small discrepancy 
between the new value of the scattering length 0.195m" 1 , with that obtained from UCPTH 
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approach 0.23m" 1 , is due to the fact that r$ given by Eq. (]2"1 ) is not much larger than 
unity. 

For a more complete 2 loop calculation of the scalar pion form factor, see the recent 
work of Gasser and Mei/3ner and Dognohue et. al. [H ,||31||. It should be remarked that 
the sum of the last 2 terms on the R.H.S. of Eq. (|3[) is not necessarily real, one must take 
their real part in the calculation. 

Eq. (|23| ) develops a pole at s g = —2.7GeV 2 which is a ghost and has to be removed 
by multiplying the R.H.S. of this Eq. a factor (1 — s/s g ) which does not affect the phase 
theorem, but does change the modulus of the form factor. Unlike the vector pion form 
factor, this ghost is near to the physical region and has to be taken into account. At 
s = m 2 K , where mx is the Kaon mass, the correction enhances the modulus of the scalar 
pion form factor by 10%. 

To end this section let us remark that in a recent study, the role of the left hand cut for 



the scalar form factor has recently been studied in details by Boglione and Pennington |32 
in connection with the criticism of the recent work of Dobado and Pelaez [|17| . They found 
that this method is an unreliable way of determining the chiral parameters but a good fit 
to the data can be done with an appropriate adjustment of the subtraction constant. We 
agree with this statement, but would like to point out that the situation for the vector 
pion form factor is different as discussed above. 



4 UCPTH and Vector Meson Dominance 

In the zero width approximation and at the tree level, the vector pion form factor can be 
written as: 

f o o 

V°{s) = J -^^ (31) 

Sp s 

where /° and g Q are, respectively, the photon-p coupling ( multiplying with e) and 
the strong pirir coupling with f^g° p7tlv /m 2 p = 1. At the tree level, these coupling constants 
should be independent of s, but loop corrections to the 2 and 3 point functions can 
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introduce the s dependence. We now want to make a finite width correction to this Eq. 
by introducing the self energy correction for the p propagator. Summing the geometric 
series for the p self energy and performing the mass and wave-function renormalisations, 
we have: 

Via) = IM9pe;M_ (32) 

S p - S- TC p (s) 



with 



p / \ (s - s p ) 2 p f°° Imn p (z) - Imir p {s p ) - (z - s p )Imn'(s p ) 

Re-n p {s) = P —P / — r dz (33) 

71 JAm\ {Z — Sp) Z {Z — S) 



and 



Amlf' 2 

where P stands for the principal part integration. The R.H.S of Eq. (BBf) could 



I™ P (s) = ^ jj? 9%«(s) (34) 



be straightforwardly expressed in terms of the function H^{s) and its derivatives, if 
g pwn (s) was a constant. In the usual VMD model, g p7T7T (s) and f p (s) are constants hence 
Qp-Kirfpl ( s p~ tt(0)) = 1. If we neglected the finite width correction, we would have = s p 
and hence the KSRF relation becomes g p7Tn = SR 1 /{2f^). Let us now use a theorem stat- 
ing that two real analytic functions, having the same discontinuity, can only differ by a 
real polynomial. By comparing Eqs. (|33|), ([34]) with Eq. ( |T2D we see that they have the 
same imaginary part and satisfy the same boundary condition and high energy behavior, 
hence they must be identical. In sum, we have the equivalence between the UCPTH and 
the VMD. 

We could improve the Vector Meson Dominance result by introducing the inelastic 
effect in the vertex correction, e.g. the p — > utt to 2tt state. This correction can be phe- 
nomenological represented by a real polynomial in s, and in order to fit the experimental 
r.m.s radius, we could choose this polynomial as (1 + 0.15s /s p ) which yields the same 
result as the calculation of the pion form factor with the inelastic contribution, Eq. (|r6|). 
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5 Generalisation to SU(3) 



It is straightforward to generalise the CPTH and UCPTH to SU(3). We have to calculate 
in addition, the KK contribution to the pion form factor. Similarly we calculate the irK 
vector form factors with the ttK and r)K contributions and also the KK and r]K form 
factors. Because the scalar i]K form factor is negligible due to the small difference of 
the K and r\ masses; the scalar ttK scalar form factor is previously well calculated, we 



refer the readers to the original calculation [IS]. Because we later want to calculate the r 



decay rate, let us now calculate the charge current (isovector) matrix element which we 
denote now by Vij, where i and j denote the pseudoscalars. The normalisation of at 



zero momentum transfer is given by the Ademollo-Gatto theorem ||33|| . Using UCPTH, 
the 7T + 7r° form factor is given by: 



WOO = V2{1 - — - ^^[{s - Aml)HUs) + ^ + \{{s-Am 2 K )H KR {s) + 

sri 96tH/^ 3 2 3 

(35) 

where the subscripts in the function H refer to the intermediate states defining this 
function. Similarly we calculate the K + K° isovector vector form factor: 



v K+K o{s) = {i - — - zjnnt&i* ~ 40n m (s) + ^) + ((* - ^ k )h kR {s) + ^f)]}- 1 

Siii 96tt - z JZ 2 3 3 

(36) 

It should be noticed at this point that in the VMD calculation, the V 7T + 7T o(s) and 
V K+K - (s) should be equal and is given by the expression for 14+vr°(s). 

The problem of Kir form factor is discussed in some details in reference We quote 
their results on the calculation of the contribution of Kn intermediate state and add to 
it the contribution of the Kr\ intermediate state (or loop). To simplify our writing, let us 
define the following function for the unequal mass case: 

Hijis, ml m)) = s 2 f ^ z*(T-T-l) d * ^ 
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where A(s,m 2 ,m 2 ) = (s — (rrii + rrij) 2 )(s — (mi — rrij) 2 ). The function Hij(s,m 2 ,m 2 ) 
can be expressed in terms of the logarithm function as was done in the references^, 
In term of this function, we can write down now the expression for the K°ir + and K + tt° 
vector form factors: 

s 1 

V K o n +{s) = {1 — (H wK (s, ml, m 2 K ) + H vK (s, m 2 m 2 K ))}- 1 (38) 

s Rs vm^f 2 

1 s 1 - 

V^ +7r o(s) = -={l — (H nK (s,ml,m 2 K ) + H vK (s,m 2 m 2 K ))}- [ (39) 

and the Kt] form factor is given by: 
[3 s 1 

v^+^(s) = - — - 1287r 2j 2 (^( g ' m ^ + ^V( s > (40) 

Strictly speaking, all above equations are derived not by the inverse amplitude method 
but rather by the Pade approximant method. It should also be mentioned that the 
K + K° and the K + 7] form factors cannot be calculated at their threshold values by CPTH 
technique because their thresholds are above the resonant masses. 

The vector pion r.m.s. radius squared is calculated to be < ry >= 0.40 ± 0.01/m 2 
compared with the experimental value < r v >= 0.439 ± .008/m 2 []30|1 . The vector Ktc 
r.m.s. radius squared is calculated to be < ry >= 0.27 ± 0.008/m 2 compared with the 
experimental value < r\ >= 0.34 ± .03/m 2 |33[. The agreement between the theory and 
the experimental data is not satisfactory (see below). 

Including the inelastic effect due to the K + K° intermediate state in the pion form 
factor calculation, changes little the numerical result. In a more elaborate calculation 



22| where the nn phase shifts were calculated first which includes the left hand cut 
contribution to the partial wave strong irn amplitudes and then compute the vector pion 
form factor using Eq. (||), there was very little change from the result of the present 
calculation. The contribution of the ujtt inelastic effect in the unitarity equation for the 
form factor is, however important [^7j . 



2(i 



The effect of the left hand cut in the P-wave elastic amplitude is, however, more 
important in the Vector Kn form factor calculation JT^, [22j . 



6 Application to Tau Decays 



Obviously CPTH cannot be used for calculations in r decays except for a very small region 
of the phase space where the one- loop calculation is valid . In this region of the phase 
space, the CPTH results are the same as those obtained by the UCPTH, but the latter 
method has a much wider range of validity. Because UCPTH can handle resonance, we 
can use this method to calculate the r hadronic decays. Let us now calculate the decay 
of the lepton r into 2 pseudoscalars i and j. Let us define the ratio of the rate of 
r — > PiPjV to that of r — > evv. Taking only into account of the vector form factor 
contribution, we have: 



- it coM : > r A(s ' m :;'" ?)3/2 (i + - i w ? * <«> 



4m r sin 2 6 c Js t s3 

where 8 C is the Cabbibo angle with sin 2 6 c = (0.222) 2 , St is the square of the threshold 
energy and M is the r lepton mass. Let us use the following notation Rij(ab) to denote 
the ratio R defined in Eq. (|4ll) with intermediate states a and b and R n K means the sum 
of the decay rate into tc + K° and tt°K + . Using Eqs. ([35l - f40| ), we have: 



R^+ n o (7T7r) = 1.05 i? 7r + 7r o(7T7r, KK) = 1.03 R-k+ttO 



exp 



1.38 ± .02 



(42) 



R k+R0 {-kix) = 0.0042 R K+R0 (irir, KK) = 0.0031 R K+R o \ e x P = -0075 ± .002 (43) 



R^ k {-kK) = 0.050 R nK (irK,r}K) = 0.048 R nK 



exp 



0.065 ± .008 



(44) 



R v k(ttK) = 5.1. 10" 4 R vK (ttK,t]K) = 4.1. 10" 4 R vK \ exp = (1.3 ± .04)10" 3 (45) 



The experimental data are taken from the recent paper of the CLEO group [34| and 
the Particle Data Group |35| . 
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It is seen that our calculation for the p and K* decays are too small by 35% compared 
with the experimental rates. This is expected because we neglect the inelastic contribution 
and essentially use the threshold parameters, the r.m.s. radii as inputs to extrapolate to 
the resonance region which is far away. 

Similarly the K + K° and K + t] decays are too low by a factor of 2. (The finite width 
correction is very small for the K*\ the KSRF relation gives IV* = 55.2 MeV and the 
finite width correction with ttK loop gives 55.6 MeV and with 7iK and i]K loops give 
53.4 MeV). The scalar form factor contribution to the above decay rates are small. Their 



maximum contribution is inR^x which amounts only to 3-4% of the calculated rate Jl9 



If we are willing to introduce more parameters in our calculation in order to fit the data 
on the top of the resonances and want only to predict the r.m.s. radii and the inelastic 
form factors K + K° and K + r], we can certainly do a better job. This can be done by using 
the polynomial in Eq. (|^) as explained above. 

Because on the top of the p resonance the peak value of the form factor is too low by 
30%, we can set P n (s) = (1 + 0.15s/s p ) in order to have the correct peak value, and hence 
we multiply the right hand side of Eqs. (|35| - ^0D by this factor. 

The vector pion r.m.s. radius squared is now < ry >= 0.46 ± 0.01/m 2 compared 
with the experimental value < r v >= 0.439 ± .008/m 2 |30|]. The vector Kir r.m.s. radius 



squared is now < r v >= 0.31±0.008/m 2 compared with the experimental value < r\ >- 
0.34 ± .03 fin 2 [j35| . There is now a good agreement with the data. 



The r decay rates are now: 



R^o(-K7r) = 1.41 R n + n o(7m,KK) = 1.38 R n+7T o \ exp = 1.38 ± .02 (46) 

Rk+koM = 0.067 R K+R0 {inT,KK) = 0.064 R K+K o U P = .075 ± .02 (47) 

R 7rK (nK) = 0.067 R^k^R, r)K) = 0.064 R nK \ exp = 0.065 ± .008 (48) 

R vK (7iK) = 8.9.10~ 4 R vK {irK,r]K) = 7.2. 10" 4 R vK \ exp = (1.3 ± .04)10~ 3 (49) 
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It is seen that the agreement with the data is much better which is not surprising 
because we do not have such a long range of energy to extrapolate, from the resonant 
energy to the r\K and KK energy available in the r decay. These results are changed 
slightly when we change the p mass to 0.773 GeV and the correction factor to (1 + 
0.13s/s p ). 

7 Conclusion 

We have presented here the study of the form factor problem using the UCPTH approach. 
The main feature of this method is that, at low energy, it coincides with that derives from 
the CPTH method. Because of the unitarity constraint, we are forced to rewrite the 
CPTH results as an infinite geometric series using either the inverse amplitude method 
or the Pade method. The former one is more general than the latter. In the simplest 
approximation to the one loop CPTH, both method yields the same result. They enable 
us to extend the analysis of CPTH to the resonant region and beyond without introducing 
more parameters. 

On the top of the p and K* masses, our calculation on the magnitudes of the form 
factors are too low and only accurate to 15% in the amplitudes. The p width is correctly 
predicted within a few percents, while that of the K* is 10 % too high compared with 
the data. The r.m.s. radii are in agreement with the data. Extending our calculation 
to above 1 GeV, then the accuracy becomes worse. The r to p and K* decay branching 
ratios are accurate to about 35% but the r to K + K° and K + rj are too low by a factor of 
2. 

In order to improve our predictions, we modify our UCPTH results which are obtained 
based on the assumption of the elastic unitarity, by incorporating the inelastic effect 
through the phenomenological polynomial ambiguity of the Omnes function. Mutiplying 
the UCPTH results on the form factors by the factor (1 + 0.15s/s p ) the leptonic widths 
of the p and K* are in agreement with the data. The prediction of the branching ratios 
r — > K + K°h / and K + rju are also in agreement with the data. The corresponding rms 
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radii are also in agreement with the data. 
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Figure Captions 

Fig.l : The square of the modulus of the vector pion form factor V(s) as a function 
of the energy squared s. The experimental data are taken from the review article by D. 



Morgan and M. Pennington ||36|| . The solid line curve is the result of the UCPTH, Eq. (|l2l) ; 
the dashed curve is the result of the standard one loop CPTH using the same parameter 
as that of UCPTH Eq. (H) ; the dotted curve is that obtained by simulating the inelastic 
ujti contribution [[^J by multiplying the UCPTH result with the factor 1 + 0.15s/s p . 

Fig. 2 :The phase of the vector pion form factor in degrees as a function of the energy 
squared s. The unitarity relation requires that it has the same phase as that of the 
strong P-wave nir scattering. The experimental data for the P-wave phase shift are taken 
from |36] . The solid curve is the calculated phase of the pion form factor using UCPTH, 
Eq. (|I2"D ; the dashed curve is calculated with the CPTH, Eq.(|^) ; the inelastic effect due 
to the ujit contribution does no affect the phase of the UCPTH form factor below the urn 
threshold. 

Fig. 3 : Test of the convergence of the perturbation expansions. The solid line represents 
the ratio ry of the one loop vector form factor as a function of the energy squared s ; 
the long dashed line represents the corresponding ratio ry for the two loop calculation; 
the short dashed line represents the ratio r$ of the scalar form factor as a function of the 
energy squared s. The convergence of the perturbation series requires that these ratios 
to be much larger than unity. 

Fig. 4 : The modulus of the pion form factor squared calculated by the UCPTH as 
compared with the sum of the first four terms of the Taylor series expansion of the 
UCPTH calculation. 

Fig. 5 : The phases of the scalar pion form factor, in degrees, as a function of the energy 
squared s. The experimental data are taken from the review article by D. Morgan and M. 
Pennington PS] . The solid line curve is the result of the UCPTH, Eq.fl2"6|) ; the dashed 
curve is the result of the standard one loop CPTH using the same parameter as that of 
UCPTH, Eq.(PD. 
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